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1 Introduction 



The use of group theory to study bound and scattering states has found wide application 
in physics 0. A common starting point is to identify a spectrum generating algebra or 
SGA for the problem of interest. This is possible when the Hamiltonian can be expressed 
in terms of generators of some algebra G. Representation theory then can be used to iden- 
tify exactly solvable limits of the theory, which can then be translated to explicit forms of 
the Hamiltonian0. Connections between bound state or scattering state problems and 
representations of compact and non-compact groups are well understood. The remaining 
category, having to do with band structure and periodic potentials, was suggested initially 
in [[3], but remained an open problem until recently. In || it was shown that dynami- 
cal symmetry techniques and representation theory can be used to solve band structure 
problems. In particular, for the Scarf Hamiltonian, band structure arises when one uses 
the complementary series of the projective representations of su(l, 1) (and so(2, 2) for the 
extended Scarf potential). Further, these dynamical symmetries allowed simple evalua- 
tion of the transfer matrix and dispersion relations. In the present study, we consider 
the band structure problem associated with the Lame equation, which is not a dynamical 
symmetry situation, but more generally that of a SGA. We discuss several realizations 
of the Lame equation and its relation to su(2) and sw(l, 1), including a discussion of 
dynamical symmetry limits, generators, dispersion relation and transfer matrix. 

While the Lame equation at first might seem to be an obscure differential equation, 
it does find surprisingly wide application in physics. For instance, it has been shown 
that if one uses the periodic potentials in super-symmetric quantum mechanics, the su- 
per partners of the Lame potentials are distinctly different from the original except for 
n = 1 which is not self-isospectral. This provides new solvable periodic potentials 0. 
The Lame equation also appears in the topics ranging from solitons to exactly-solvable 
models. This includes associations with solutions of the periodic KdV equation^, BPS 
monopoles0, sphaleron solutions of the (l+l)-dimensional abelian Higgs model||, sine- 
Gordon solitons P], as well as relations to Calogero-Moser systems |TD[. The analysis of 



this equation has many group theoretical aspects as well. For instance, the Lame equa- 
tion can be written in terms of the composition of two first-order matrix operators, where 
the coefficients of which satisfy so(3) Nahm's equationjll]. It also arises naturally in 
the context of su(2) when one tries to separate Laplace's equation in certain coordinate 
systems fT2|, |T3f as well as in the general classification of Lie algebraic potentials |TJ[] and in 



quasi-exactly solvable sl{2) models^]. The discussion of band structure in the context 
of group theory is more recent however 0, although the system was not solved and the 
transfer matrix was not determined. 

The Jacobian form of the Lame equation is 

- + K 2 i{i + l)sn 2 x = S. (1) 
ax 2 

We will view this as a Schrodinger equation with mass M =1/2 when x is real valued. This 
equation was first studied group theoretically in [|13|]. The elliptic functions sna =sn(a|/t), 
cna = cn(a|/«), and dna = dn(a|ft) are doubly-periodic functions in the complex plane, 
of modulus k, where < k < 1. We will omit the modulus except in functions where 
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it differs from k. The complementary modulus is defined as k' — (1 — k 2 ) 1 / 2 . The 
periods of the Jacobi elliptic functions are related to the complete elliptic integrals K = 
(tt/2)F(1/2, 1/2, 1; k 2 ) and K = (tt/2)F(1/2, 1/2, 1; k' 2 ), where F is the hypergeometric 
function, by: 

sn (a + t) = sn a t = 2iK', AK + UK', AK (2) 
cn(« + r)=cna r = 4iiT, 2K + 2iK, AK 
dn (a + r) = dn a r = AiK, AK + AiK, 2K 

In the limit k = 1, the real period related to K becomes infinite, and the functions are no 
longer periodic on the real axis. Along the real axis, the potential sn 2 x is periodic and 
bounded, while along the imaginary axis, it has periodic singularities of the type 1/x 2 . 
Examples are shown in Fig. 1 for k 2 = 1/2 for potentials V(x) = sn 2 (ax + 6), with several 
choices of a and b. 

The Lame equation has also been related to the band structure problem associated 
with the Poschl- Teller potential, V p t(x) ~ 1/ cosh 2 x \TE, M, by using the expansions for the 



elliptic functions in terms of trigonometric functions found in the exercises of Whittaker 
and Watson |p77f1 . Specifically, one can start with the 1-d band structure problem 

d 2 00 a 
dx 2 fc £^ jo cosh 2 ((x - ka)/x ) 



If we write the coupling constant as 



x. 



^ (4) 







we can make the coordinate transformation 



which transforms the Hamiltonian to the Lame equation. In this article, we first discuss 
the relation of su(2) Hamiltonians to the Lame equation, followed by a discussion of band 
edges. Next we consider su(l, 1) and the relation to scattering states in the dynamical 
symmetry limits. Finally we derive the transfer matrix and dispersion relation for the 
Lame equation followed by the relation to the group theoretical SGA Hamiltonians. 

2 su(2) Realizations of the Lame Equation 
2.1 Sphero-conal Coordinates for su(2) 

In order to realize the Lame equation from su(2) or su(l, 1), we will use coordinate systems 
defined in terms of Jacobi elliptic functions. These are often refered to as sphero-conal [ jjf 



or conical|12] coordinates. One can find variations that have been considered in the past, 
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and not all provide a good description for all values of the modulus k. Consider first the 
mapping: 

k 1 

x = k sn a sn (3, y = i — cn a cn (3, z = — dn a dn (3, (6) 

k' k' 

The ranges of the angles a, (3 are defined in terms of the elliptic integrals K and K', so 
that they depend on the value of k. Specifically — 2K < a < 2K and K < (3 < K + 2%K' 
where k' 2 = 1 — k 2 , and x 2 + y 2 + z 2 = 1. The contours of constant a and (3 cover the 
sphere as illustrated in Fig. 2(a). For k — 1/2, four typical contours (for a = ±K/2 and 
(3 = K + iK'/2, K + 3iK'/2) are shown in Fig. 2(b). This parametrization can be used 
to construct the generators of su(2). A direct computation of the generators yields the 
realization (I) given in Table 1. These generators satisfy the commutation relations in the 
form [L a , L b ] = e a b c L c with Casimir invariant 



2 , r 2 , r 2 

/ d 2 d 2 



C 2 = L 2 X + L 2 + L Z 



k 2 (sii 2 a — sn 2 (3) \da 2 d(3 2 



(7) 



The su{2) algebras we discuss here have Casimir invariants which have expectation value 
— £(£ + 1). We now consider which types of quadratic Hamiltonians of the form 

ir = 5>Lj (8) 

3 

can be constructed which result in Schrodinger equations in a and (3. Specifically, when 
the coordinate representation of H contains differential opertors only of the type d 2 /da 2 
and d 2 /d(3 2 , one can use H together with C 2 to perform a separation of variables resulting 
in two independent Lame equations. We find the three forms in the top right of Table 1, 
which are simply related through the Casimir operator. To demonstrate the separation, 
consider H 3 , 



Hs = L 2 z + n 2 L 2 y 
1 



(sn 2 (3 — sn 2 a) 




(9) 



and require it to be a constant of the motion with eigenvalue H 3 = £. Using basis states 
denoted by the direct product, ^(a, (3) = ip(a)(f)((3), we arrive at two identical Lame 
equations with identical eigenvalue: 



-^ + K 2 £(£+l)sn 2 a 
da 2 

d 2 



dp' 



+ k 2 £(£ + 1) sn 2 /3 



4>(a) = Eip(a), (10) 
m = £<t>{(3)- 



Thus solution of the eigenvalue problem for H 3 = L 2 + k 2 L 2 yields the single valued 
solutions to these equations, which will correspond to the band edges. One should keep 
in mind that while a is defined on the real axis, (3 is complex. Also, while the Hamil- 
tonians H k are well defined in their algebraic form, the generators and coordinates have 
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singularities in both k — and k — 1 limits. Hence we consider a different realization 
below. 

To reduce these Hamiltonians to more conventional Schrodinger equations in a and 
f3, we would like the coordinates to be defined along the real axis. To do so, we make a 
shift (i) (3 — > (3 + K + iK followed by (ii) f3 — > —ij3. This results in the transformations: 



sn 



I") 



-dn(/3|/c'), 



-cn(/V), 



IK 



cn(/5|/c) ■ 
dn(/3|«) -> «'sn(/3|K / )- 
The new coordinate system is no longer singular, and is parameterized as 
x = sna da(P\n'), y = cna cn(/3|/c'), z = dn a sn(/3|/c'), 



(11) 



(12) 



where -2AT < a < 2AT, -AT < /3 < AT are now both real valued. Contours of constant 
a and /3 are shown in Fig. 2(c) for a = ±K/2 and /3 = ±K'/2. If we construct the 
generators of su(2) in this basis, we obtain the realization (II) given in Table 1. The 
Casimir invariant now has the form 



C 2 = Ll + Ll + Ll 



-1 



d 2 d 2 
+ 



K 2 cn 2 a + K ,2 cn 2 (f3\K') \da 2 d(3 



The quadratic Hamiltonians which give rise to Schrodinger equations when combined with 
C*2 are given in the lower right of Table 1. Consider again H 3 = L 2 + k 2 L 2 . If we set 
H 3 = £ and C 2 = —i{t + 1), we obtain 



d 2 



da 2 



+ K 2 i(£ + l)sn 2 a 



dp 2 



+ K ,2 £(£ + l)sn 2 ((3\ K ') 



if} (a) = £ip(a), 



(13) 
(14) 



In this case the coordinates are both real valued. For integer £, and both potentials are 
non-negative, we see that < £ < £(£+ 1). Further the ground state of one Hamiltonian 
corresponds the the highest energy state of the other, and vice-versa. (Similar equations 
are obtained if we use H 2 or Hi.) Before we discuss some of the general properties of these 
equations, we consider first the dynamical symmetry limits which are when k — or 1. 



2.2 k = Limit 

In the limit k — 0, and H 3 = L 2 = m 2 , so that £ = m 2 , and the coordinate Hamiltonians 
are: 

d 2 



da 2 

d 2 _ 1(1+1) 
'dp 2 ~ cosh 2 f3 



m = -em- 



(15) 
(16) 



5 



The generators in this limit have the simpler form: 



h 



=F cosh (3 7^- + % sinh f3-^- 
op oa 



.d_ 

da 



(17) 



The first equation is free motion, while the second corresponds to bound states of the 
Poschl- Teller potential. For the discrete representations of su(2), one obtains the bound 
state spectrum —S = —m 2 . The wavefunctions are of the form: 



^(a, (3; k = 0) ~ PP(tanh f3) exp(±ima) 
with —Ti < a < Ti, — oo < (3 < oo. 



m = — t 



0,1,2, 



'11 



2.3 k = 1 Limit 

In the limit k = 1, H 3 = L 2 — L 2 X1 the equations reduce to 
d 2 £{£ + lY 



da 2 



cosh a 



d 2 



dp- 



if>(a) = [S -£{£ + l))ip{a) = S'tp{a) 

4{p) = W + 1) - s\m = -£W) 



(19) 
(20) 



The generators are the same as in the k = limit, with a and f3 interchanged. The second 
equation is now free motion, while the first corresponds to bound states of the Poschl- 
Teller potential with a shifted eigenvalue. For the discrete representations of su(2), one 
obtains the bound state spectrum £' = —m 2 . The wavefunctions are of the form: 



^(a,(3;K = 1) ~ P £ m (tanha) exp(=b"m/3), 
with — oo < a < oo, — n < (3 < n. 



m = —I 



0,1,2, 



(21) 



2.4 Band Edges and su(2) 

For general values of k, the Lame Hamiltonians are periodic (as seen in Fig. 1), and will 
have band structure. As the discrete representations of su{2) correspond to single-valued 
wavefunctions, the discrete representations can at most describe the band edges. Since 
we have an algebraic realization of the Hamiltonians whose spectrum corresponds to the 
Lame equations, we can obtain the eigenvalues of the Hamiltonian Hk (and hence the 
band edges) by a direct diagonalization in the spherical harmonic basis \£m), which can 
then be related to the results for the elliptic basis through the coordinate transformations. 
The resulting functions are doubly periodic solutions of Lame's equation known as Lame 
polynomials ||18||. These polynomials are of the form sn a xcn b xdifxF p (sn 2 x), where a,b,c = 
0, 1 and a + b + c + 2p = £. Here F p (z) is a polynomial in z of order p. A discussion of 



6 



these functions can be found in jyj [u], [LJ|. For I — 1, there are three eigenstates of H 3 
given by 

= |10) ~ snadn(/5|K') 
#2(0, #«) = -J=[|ll) + |l-l>] ~ cnacn(/V) (22) 

tf 3 (a,/3;«) = dnasn(/V) 

with eigenvalues i?i = 1 + k 2 , E 2 = 1, E3 = k 2 . In the dynamical symmetry limits, these 
three states reduce to those discussed above. 
For £ — 2, we have 

tfi(a,/?;«;) = -p[|22) - |2-2)] ~cnadnasn(^|K')cn(/3|«') 



V2 

ty 2 (a,(3;K) = -=[121) + |2- 1)] ~snacnacn( / 9|«0dn( / 9|« / ) 
v2 

* 3 (q;,/3;k) = -p [|21) - |2 - 1)] ~ snadnasn(/V)<M/V) (23) 
v2 



'2 2/+(«) - 1 - k 2 

. 3 

;i-/ + ( K )sn 2 a)(l-/_(^)sn 2 (/?|«:')) 
^2 2/_(/t) - 1 - 

3 1-K 2 



¥ 4 (a,/3;«) - |22) + |2-2) + t/- J+ 7 2 1 20> 



* 5 (a,/3;K) ~ |22) + |2-2) + J- ^/ ^ 120> 



(l-/_( K )sn 2 a)(l-/ + ( K ')sn 2 (/3|«')) 



with eigenvalues 



E x = 1 + K 2 , E 2 = 4 + k 2 , E 3 = 1+4k 2 , (24) 
£4 = 2/ + (k), £5 = 2/_(«), (25) 



where /±(k) = ^1 + k 2 ± vT— k 2 k' 2 ) . In general there will be a relation between the 
eigenstates |£m) of the 2£ + 1 band edges and products of Lame polynomials. The band 
edges for £=1,2 are shown in Fig. 3 as a function of k 2 . The bands are indicated by the 
shaded regions. The dashed line indicates the height of the potential. One can see that 
the bands merge at k = 0, and pass to the bound and scattering states of the Poschl- Teller 
potential for k = 1. 



3 su(l,l) Realizations of the Lame Equation 
3.1 Elliptic Parametrization of su(l,l) 

In order to discuss the eigenstates in the band, and develop the dispersion relation, we 
require more general representations. In the spirit of recent work on the Scarf potential 
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where it was shown that one can use a su(l, 1) dynamical symmetry to analytically solve 
for the dispersion relation^], states and transfer matrix, we consider transforming our 
generators to su(l, 1). Consider the transformation of our previous coordinates associated 
with x — > —ix, y — > —iy : 



x = —in sn a sn (3, y = — cn a cn (3, z = — dn a dn /3, 

K 1 



K 



(26) 



where < a < 4K, < (3 < iK |TJ. This now parametrizes the hyperbolic surface 
z 2 — x 2 — y 2 = 1 illustrated in Fig. 4(a) for k 2 = 1/2. Contours are shown for selected 
values of a and f3 in Fig. 4(b). The generators of su(l, 1) algebra are given in Table 2 as 
realization (I), and satisfy the commutation relations [L x ,L y ] = —L z , [L y ,L z ] = L x , and 
[L z , L x ] = L y . The Casimir invariant now has the form 



Co 



d 2 



d 2 



/t 2 (sn 2 a 



sn 2 (3) \da 2 dp 2 



(27) 



As in the su(2) case, there are three forms of bilinear Hamiltonians which lead to Schrodinger 
equations in a and (3. These are given in the top right of Table 2. If we choose 
H 3 = L 2 — k 2 L 2 x = S, together with the Casimir invariant Ci = —£(£ +1), we obtain 
the decoupled Lame equations: 



d 2 



da 2 
'd(3 2 



+ k 2 £(£ + l)sn 2 a 



+ k 2 £(£ + l)sn 2 /? 



-em- 



(28) 
(29) 



Again similar results are obtained if we use Hi or H 2 , the only difference arising in the 
definition of the eigenvalue. Note that the range of j3 is along the imaginary axis. This 
realization is problematic in the k = and k = 1 limits since the generators become 
singular. Consequently we consider a slightly different realization of su(l, 1) below. 

3.2 Another su(l,l) Realization 

While we do not have a coordinate system which avoids the singularities at K — 1, it is 
possible to at least allow a study of the scattering states when k = 0. To do this we make 
the transformations (i) (3 — > (3 + K + iK followed by (ii) (3 — > —i(3. As a result: 



sn{p\K) 

cn.((3\n) 
dn((3\K) 



-dn(/V) 



K 

— cn( 

IK 

ft'snf/ 



k 



k 



(30) 



This results in coordinates 13 



x = — i sn a dn(/5|/t'), y = —icna cn( 



dna sn( 3\k' 



(31) 



8 



where < a < AK, —iK < (3 < —iK + K' , also satisfying z 2 — x 2 — y 2 = I. Typical 
contours for this parametrization are shown in Fig. 4(c). The generators are given in 
Table 2 as realization (II) and satisfy [L x ,L y ] = —L z , [L y ,L z ] = L x , and [L Z ,L X ] = L y 
with Casimir invariant 



C = L\-L 2 X 



L 2 = 



d 2 d 2 
+ 



K 2 cn 2 a + K /2 cn 2 (f3\K') \da 2 d(3 2 



(32) 



The Hamiltonians which result in seperable Hamiltonians are given in the bottom 
right of Table 2. Using K 1 = L 2 X + K ,2 L 2 y = S and C 2 = -1(1 + 1) leads to the two Lame 
equations: 



d 2 



da 2 



+ K 2 i(£ + l)sn 2 a 



d 2 



dp' 



+ k' 2 £(£ + 1)sii 2 {P\k' 



(33) 



It should be kept in mind that (3 is still complex valued. Never the less this realization 
allows an analysis of the k = dynamical symmetry. 

3.3 k = Limit 

This realization is not singular in the k = limit (although it is in the k — 1 case). 
Taking k — 0, and shifting (3 to be on the real axis by (3 = 9 — in/2, we find 



L± 

L, 



( d d N 

±e ±ia =F sinh 9— + i cosh 9— 
\ o9 oa / 

d_ 

da 



(34) 



To recover the usual su(l, 1) commutation relations [I Z) I±] = ±I±, [I+,I-\ = —2I Z , we 
then identify I± = rkL± and I z = %L Z . It is convient to perform the transformations: 
9^9 — in/2, tanh# — > cos#, followed by a similarity transformation f(9) = \J sin 6, and 
9 — > i9. Then we have the form: 



da~ 2 ^ a) 



d 2 m 2 -l/4 
~d9 2+ sinh 2 # 



m 2 ^jj(a), 



(35) 



The principal series I = —1/2+ip of the projective representations of su(l, 1) now describe 
these scattering states, and the eigenfunctions are of the form: 



~ yji sinh 9PP (cosh 9)e ±l 



±ima 



m£8, £= -1/2 + ip, p>0. 



(36) 
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4 Band Structure of the Lame Hamiltonian 



We now focus on the properties of the Lame equation in the form of Eq. (1). As we 
are interested in band structure, the eigenstates must satisfy Bloch's theorem. For a 
potential which is periodic with period a, V(x + a) = V(x), the wavefunctions must be 
of the form 

^k(x) = v>k(x) exp[— ikx], (37) 

where k is the wavenumber and Uk(x) has the periodicity of the lattice: Uk(x + a) = Uk(x). 
As the eigenstates ^(x) are not periodic, the doubly periodic solutions of the Lame 
equation do not play a role for energies in the band. Rather, we look to the more general 



class of solutions expressed in terms of Jacobi theta functions ||17||. Starting with the 
Hamiltonian 

" d 2 

+ K 2 £(£ + l)sn 2 (x\K) 



Hip 



dx 2 



ip(x) = £ip(x) 



(3? 



the solutions for positive integer I are given parametrically by 



i>(x) = n 



71=1 



H(x + a 



n}_ -xZ(atn) 



6(x) 



(39) 



where Ti and 9 are theta functions, and a\, (%2, ■■■olh are constants determined by the 
constraints: 



£ 




E 

n=l 

e 

E 

P =i 



i 



sn 2 a r 



cna n dna n /sna r , 



.71=1 



sna p cna p diaap + sna„cna n dna r> 
sn 2 a p — sn 2 a n 



(p ^ n) 



(40) 
(41) 



If this solution is not doubly periodic, a second solution is 



if) (x) 



n 

n=l 



H(x a n ) x z(a n ) 
9(x) 



(42) 



We can then identify the dispersion relation by putting these wavefunctions in Bloch form 
and by using the periodicity of the theta functions to extract Uk(x). We find 



k{£) 



J2Z(a n \K 2 ) + 



71=1 



2K' 



(43) 



We will start with the case of I = 1, where simple anlytic results can be obtained. We 
then derive the transfer matrix for the general case of integer £. 



4.1 t = 1 Results 

Since there is only one parameter a, the constraint equation is simply 

dn 2 « = £ — k 2 . 



(44) 
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The condition that the dispersion relation is real, ReZ(a\K 2 ) = 0, results in two energy 
bands given by 

k 2 <£<1, 1 + k 2 <£. (45) 

These are shown in Fig. 3 (top). In the lower band, a has the form a — K + irj, where rj 
ranges from K' at £ — k, 2 , to at £ — 1. In the upper band, a = ir], where rj ranges from 
at £ = 1 + k 2 , to K' as £ — > oo. This path traced out by the parameter a as a function 
energy £ is shown schematically in Fig. 5. The upper and lower sides correspond to band 
gaps while the right and left edges are the energy bands. Using the specific forms of a for 
each band, the dispersion relation becomes 

k(s) = I -^K 2 ) + £(i - + ^r^ 1 - 2 < z < i. m 

This is plotted in Fig. 6 (a) for the case of I = 1. The momentum k is plotted up to the 
edge of the Brillouin zone, which is k — n/2K. In the figure we use k 2 = 1/2, so that 
the band edges are £ = 1/2, 1 and 3/2. (The analogous behavior for £ = 2 is indicated in 
Fig. 6(b).) The solution of the Lame equation in Bloch form is now 



Hl ^ v) exp(i7rx/2K)} exp(-ikx) k 2 < £ < 1, 



^ggM ^ V (mx/2K)] exp(-ikx) 1 + k 2 < £ 



(47) 



The component of the wavefunction in square brackets can be checked to be periodic with 
the periodicity of the direct lattice: x — > x + 2K. 

The dispersion relation displays the desired limits. One can see that k(£) — > as 
£ — > k 2 , and k{£) — > n/2K as £ — > 1. Further, as the modulus of the elliptic function 
vanishes, k 2 — > 0, the Hamiltonian becomes that of a free system, and we find fc(£) — > 
as desired. In the Poschl- Teller limit, k 2 — > 1 , the band vanishes, and we find — > 0. 
Similar results hold for the upper band as well. 

From the dispersion relation we can compute the group velocity and effective mass. 
To do so, we use the relation between the zeta function and the elliptic integrals of the 
first and second kind 

Z(a) = E(a) - (48) 
K 

where E{k 2 ) = f F(-l/2, 1/2; 1; k 2 ) and K = f F(l/2, 1/2; 1; k 2 ) are the complete elliptic 
integrals and E(a) is incomplete. Then, the group velocity is given by: 

i _- d m_ (49 ) 

v~ d£ 1 yj 

or 

2 X /(1 -£)(£- K 2 )(l + K 2 -£) 

We plot v 2 as a function of energy in Fig. 7 for several values of k. As k approaches 
zero, the energy gaps vanish, and the group velocity approaches the free particle limit 
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E = Mu 2 /2 = v 2 / 4 (dot-dashed line). As k approaches unity, the lower band vanishes 
becoming a bound state, and the group velocity is only non-vanishing for the continuum 
states of the Poschl- Teller potential with £ > 2. 
The effective mass M* is determined by 



1 dv 
~M* ~ dk 

,(£ - « 2 )(1 + k 2 -£) + (£- 1)(1 + k 2 -£)-(£- 1)(£ - k 2 ) 



= -2- 



(£-m-K 2 ){i + K 2 -£) h 

+ 4 77. 2 W^hi (5 ^ 

[£-K 2 - -^f 

We plot 1/M* in Fig. 8 for selected values of k. One can see that as k 2 — > 0, the gaps 
vanish and M* — > M = 1/2 as expected for the free particle. 

4.2 Transfer Matrix for the Lame Hamiltonian 

The general form of the transfer matrix is computed using the definitions in the Appendix. 
Using the wavefunctions and Eq. (Al), we have 

^ sna n dna n .£ir 

r = ~g> K) - cna n ] ~ l Yk- (52) 



Consequently, we see that 



r + i k(£) = j2 Snandnan - (53) 



n=l CnQ r 



The transfer matrix then has the form 
T = 



cos 2k(£)K i (ELi ^^g 2 *) _1 ^ 2A;(£)^ \ m 

»(En=i sn "^° a " )sin2A;(g)ir cos2A;(^)ir /' 



(Note that this is in the form of Eq. (A. 3) rather than (A. 5).) While this expression 
requires knowledge of the parameters a n , one can obtain various limits of this for the free 
particle and Poschl- Teller potentials. 



4.3 The k = Free Particle Limit 

We start first with the I = 1 case. Taking k = in our transfer matrix, we obtain for the 
upper and lower bands 



cosh2k(£)K i^— smh2k(£)K \ 



_ V ^ snadna V ^ fee:") 

• - | ,- snadna sinh2A; ^^ CO sll 2fc(£) J K" )' { ' 
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According to Eq. (pEOD, 



sn a dn a 



cn a 



Taking k = 0, we have r — >■ 0, so that 



;i + /t 2 -g)(g-/t 2 
(f-i) 



.snadna /— 

= v£. 

cn a 



(56) 



(57) 



The transfer matrix becomes that for a free particle, given by 

T= I 



COS 71 sm7r 
£ 1//2 sin ixyfS cosTX\fE 



For the general case of integer £, we start with the constraint equations 
limit K — > 0. We first must show that 



r + ) = ^ tana n — > iV~£. 



(58) 

and take the 
(59) 



71=1 



Then, from the definition of r in the appendix, it must vanish for free motion, so that 
r — > implies ) — > Then we recover the free particle transfer matrix. While we 
can show that the sum in Eq. (|59|) tends to i\f£ for small £ on a case by case basis, we 
do not yet have a general proof. However, since the transfer matrix must be that of a free 
particle in this limit, it is clear that Eq. (|5Qf) must hold, and we can use this instead to 
provide an additional relation among the parameters a n . 



4.4 The k = 1 Poschl- Teller Limit 



In the limit k = 1, the Hamiltonian H becomes the Poschl- Teller Hamiltonian, and our 
transfer matrix should reduce to that case. We start first with £ = 1 and examine the 
upper band. (The lower band becomes degenerate at K — 1). In this limit, K oo and 
K' — > 7r/2. For the upper band, 

m ^^Zr2 + ——, (60) 



From Eq. (p6|), we have in the k = 1 limit: 

sn a dn a , r- , . 

= tanh a = - 2. 61 

cn a 

The asymptotic form of the transfer matrix (as K — > oo) becomes: 

/ cos(2ATy/g^2 + 2m ) (£ - 2)' 1 / 2 sin (2Ky^ ^2 + 2zq) \ , , 

I - (£-2) 1 / 2 sin(2ir v / ^ r 2 + 2ia) cos{2K^f£^2 + 2za) J" ^ j 
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We must now change this form of the transfer matrix, defined for periodic potentials T, 
to the form used in the Poschl- Teller case which is not periodic, T (see Appendix). Using 
k(S) = y/E=2 , we find: 



T 



K=l 



T{ik)r(l+ik) 
F(2+ik)T(ik-l) 







r(-ik)r(i-ik) 
r(2-ifc)r(-ifc-i) 



(63) 



For the case of arbitrary £, we have the more general relations: 

i 



E - E tanha n = iy/£ -£(£+!). 

n=l ana n n= i 



(64) 



The dispersion relation becomes 

i 

k(£) — > —i (tanh a n 



71=1 



K J+ 2K 



2fT a 



£-£(£+!) + £— + z^E««- 



(65) 



71=1 



To compute the form of the transfer matrix for the Poschl- Teller potential, we use the 
notation 



T 



F G 

G* F* 



where, from the appendix, we have: 



F = (-) £ exp 



-2E 



a,. 



71=1 



n 



tanha„ — 1 



~ x tanh a n + 1 

This can be further simplified using the limiting form of the constraint equations 



(66) 



(67) 



S 




E 



1 — tanh a r 



Ln=l 

e 

E 

p=l(p^n) 



tanh a r 



+ £ coth 2 



ot,. 



71=1 



1 — tanh 2 a n ) tanha n + (1 — tanh 2 a p ) tanha f 
tanh 2 a p — tanh 2 a n 



(6* 



For small £ (£ = 1, 2) we can solve these equations to show that 

r(ifc)r(i + ik) 



Y{l+ik + £)Y{ik 



G = 0. 



(69) 



This reduces to the correct form of the transfer matrix. For arbitrary £, we have not been 
able to solve these equations, but the result must still hold, since it is governed by the 
form of the Hamiltonian. Again, we can use these relations to provide additional relations 
between the parameters a n . 
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5 Conclusions 



We have presented a group theoretical analysis of the Lame equation, which is an example 
of a SGA band structure problem for su(2) and su(l, 1). We have computed the disper- 
sion relation and transfer matrix, and discussed the limiting dynamical symmetry limits 
of these results which correspond to the Poschl- Teller and free particle Hamiltonians. Be- 
cause the general Hamiltonian is not of the dynamical symmetry the spectrum cannot 
be obtained in closed form. Never the less, a diagonalization of Hamiltonians which are 
bilinear in the angular momentum generators will provide the general solution. There are 
still many open questions associtated with the group theoretical treatment of the Lame 
equation. It would be nice to develop a su(l, 1) parametrization which is non-singular for 
all values of k. In addition, the Scarf and Mathieu equation limits would be interesting 
to realize more expliticly in the transfer matrix and dispersion relations. Finally, the 
diagonalization of the algebraic Hamiltonians in the continuum su(l, 1) bases would be 
interesting to study. It is clear that the result of the diagonalization must yield the same 
transcendental equations for the parameters a n , but their origin would be different. 



Appendix: Form of the transfer matrix for periodic 
potentials 

For symmetric periodic potentials with period r, we can derive a general formular of the 
transfer matrix. Suppose for a specific energy E, we have two bloch solutions u + (x)e lkx , 
U-(x)e~ tkx , where x is the coordinate, k = k(S) is the dispersion relation. Since the 
potential is symmetric, we can define pl|| 

(A.2) 



«+(r/2) _ u'_(r/2) 



u + (r/2) w_(r/2)' 

Then the transfer matrix is 

_ / cos kr ttt— sinfcr \ , . 

T = \ .... v . , lk+r , . (A.3) 
V i(ik + r) sin Act cos/ct / 

The forms of transfer matrices used for periodic and non-periodic potentials is different. 
In the discussion of the Lame equation, the k — > 1 limit takes a periodic potential to a 
non-periodic one, so that we require the transformations that take us from one standard 
form to the other. If we express the transfer matrix T for the periodic potential as 

/ Re(Fexp(z£;T) + G) ±Im(F exp(ifcr) + G) \ , . 

{ klm(F exp(ifcr) - G) Re(F exp(ikr) — G) )' 1 } 



then that of the non-periodic limit will have the form 

F G 

G* F* 



r = ( Z t ) ■ (A.5) 
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In this notation, k = k{£) is the dispersion relation, and r is the period of the periodic 
potential which tends to oo. 
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Figure Captions 



Figure 1. Various forms of the Lame potential which are real valued. We plot V(x) for modulus 
K 2 = 1/2 given by (i) sn 2 (x|/t) (solid), (it) sn 2 (x + iK'\k) (dots), (Hi) sn 2 (x + K + 
iK'\k) (dashes), (iv) sn 2 (ix\K) (dot-dashes). For other values of x, the potential is 
periodic, but generally complex. 

Figure 2. Coordinate systems for su(2). (a) Parametrization of the sphere, (b) Contours of 
constant a and (3. (c) Same as (b) but for the second parametrization. 

Figure 3. Band edges from the su(2) realization of the Lame equation as a function of k 2 for 
£ = 1 (top) and £ = 2 (bottom). The eigenvalues are indicated. The bands are given 
by the shaded regions. The dashed line indicates the value of the energy which is 
at the maximum of the potential V(x) = k 2 £(£ + l)sn 2 x 

Figure 4. Coordinate systems for su(l, 1). (a) Parametrization of the hyperboloid. (b) Con- 
tours of constant a and (5. (c) Same as (b) but for the second parametrization. 

Figure 5. Evolution of the parameter a for the £ = 1 Lame equation. The lower energy band 
has parameter a = K + irj where 77 ranges from K' at the lower end to at the 
upper end. The valence band starts with a = and grows to a = iK' as the energy 
E -> 00. 

Figure 6. (a) Dispersion relation S(k) for the £ = 1 Lame equation with k 2 = 1/2. The 
momentum is plotted in units of n/2K. (b) Analogous behavior for the £ = 2 case. 

Figure 7. Group velocity for the I — 1 Lame equation for several values of the modulus k. 

We show v 2 as a function of energy S. When k — 0, the Hamiltonian is that of a 
free particle so £ oc v 2 (dot-dashes). As k increases to unity, the lower energy band 
shrinks to a single bound state of multiplity 2 corresponding to the Poschl-Teller 
potential, and the group velocity is only non-zero in the continuum, £ > 2. 

Figure 8. Behavior of l/M* as a function of energy £ for £ = 1 Lame equation, where M* 
is the effective mass. As k — > 0, the gap closes, and M* — > M = 1/2 for this 
Hamiltonian. 
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Table 1. Realizations of SU(2) with [Li,Lj] = e^Lk and corresponding Hamiltonians which lead to the Lame equation. 



SU{2) Realization 


Hamiltonians 


(I): 

L x = i [/t(sn 2 o; — sn 2 /?)] 1 

x [sn a cn (3 dn (3d a — cnadnasn f3dp] 
L y = [Kn'(sn 2 a — sn 2 /?)] -1 

x [cn a sn (3 dn (3d a — sn a dn a cn f3dp] 
L z = i [/t'(sn 2 a: — sn 2 /?)] -1 

x [dn a sn (3 cn (3d a — sn a cn a dn f3dp] 

C 2 =I? X + I? V + I? Z = [n 2 (sn 2 p - Ba'a)]- 1 (<9 2 - <9 2 ) 


Hi =Ll + n' 2 L 2 y 

= [/t 2 (sn 2 a - sn 2 /?)] -1 (dn 2 (3d 2 a - dn 2 adfj 
H 2 = n 2 L 2 x -n' 2 L 2 z 

= [sn 2 (3 — sn 2 a>] 1 (cn 2 /?<9 2 — cn 2 adfj 
H 3 =L 2 z + K 2 L 2 y 

= [sn 2 (3 — sn 2 a] -1 (sn 2 /?<9 2 — sn 2 adfj 


(II): 

L x =[K 2 cn 2 a + K ,2 cn 2 ((3\K , )r 1 

x [k /2 sn asn(/? |/c')cn(/3 |/c')9 Q + cnadna:dn(/? /t')<9g] 

L y = [K 2 cn 2 a + K' 2 cn 2 (^|K')] _1 

x [— cnasn(/? K / )dn( / 9 |/c')9 Q + sna dna cn(/? 

L z =[« 2 cn 2 a + K' 2 cn 2 (/3|K / )]" 1 

x [— dna cn(f3\K,')dn((3\K')d a — k 2 sn a cn a sn((3 \n')dp\ 

C2 =Ll + Ll + Ll 

= [K 2 cn 2 a + K ,2 cn 2 (P\K')Y l (<9 2 + d 2 ) 


Hi =LI + k»L* 

= [n 2 cn 2 a + k' 2 cyi 2 (P\k')Y 1 (k%yi 2 (PW)d 2 a + dn 2 a<9|) 
H 2 =k 2 L 2 x -k' 2 L 2 z 

= - [K 2 cn 2 a + K^cn 2 \PW)Y l (/t' 2 cn 2 (/3|<)d 2 - /t 2 cn 2 a<9 2 ) 

= [/t 2 cn 2 a + /c^cn 2 ^')] -1 (dn 2 (/3K)<9 2 + /t 2 sn 2 a«9 2 ) 



Table 2. Realizations of SU(1, 1) and corresponding Hamiltonians which lead to the Lame equation. 



SU (1,1) Realization 


Hamiltonians 


(I): 














L x 


= i (kk' (sn 2 (3 — sn 2 a)) 1 






= /t 2 L 2 + K ' 2 L 2 








x [cn a sn (3 dn j3d a — snadnacn fldp] 






= [sn 2 a — sn 2 /3] 1 ( cn 2 /3<9 2 — cn 2 a<9|] 




Ly 


= (/t(sn z p — sn z aj) 




#2 




'(dn 2 /^ 2 - dn 2 a<9 2 ) 






x [sn a cn (3 dn (3d a — cn a dn a sn (3dp\ 






= [/t 2 (sn 2 a - sn 2 /3)}~ 




z 


= i [K , {sn 2 f3 — sn 2 a:)] 1 




o 


= Ll~ n 2 L 2 x 






x [dn (x sn 3 cn /3cL — snacnadn ddn] 






= — [sn 2 a — sn 2 (3} 


( sn 2 /3<9 2 — sn 2 a<9|) 




c 2 


= L 2 — L 2 — L 2 = — [ K 2 (sn 2 a - sn 2 /3)] -1 

,0 it/ y L \ 1 / J 


(dl - df) 

\ a p J 










(II): 














L x 


= t[K 2 cn 2 a + K ,2 cn 2 ((3\K')r 1 




#1 


= Ll + K' 2 L 2 y 








x [— k /2 sn a sn(/3 |/c')cn(/3 K/)9 a — cnadn 


adn(/V)fy] 




= - [/t 2 cn 2 a + /t /2 cn 2 


[PW)}- 1 (V 2 sn 2 (/« d\ 


+ dn 2 «(9^ 


Ly 


= *[/€ 2 cn 2 a + ^W^K)]^ 




#2 


= - LI 








x [cnasn(/3|«; / )dn(/3|«; / )9 a — snadnacn( 






= — [/t 2 cn 2 a + /t /2 cn 2 


'An')}- 1 Un 2 {(3\K')dl + K 2 sn 2 ad 2 ) 




= [K 2 ai 2 a + K Q CQ 2 ( i 9|K')]" 1 




#3 


= k 2 L 2 x + k' 2 ! 2 








x [— dna cn(/9 K')dn(/5 |/c')9 a — k 2 sn a cn a sn(/3 




= [K 2 cn 2 a + k/W^IO]" 1 (V 2 cn 2 (p>') <9 2 - 


/t 2 cn 2 o; djf) 




- /~ 2 /~ 2 /~ 2 














= [/t 2 cn 2 a + «' 2 cn 2 (/3|/c')]" 1 (<9 2 + df) 













This figure "lame-fig2.jpg" is available in "jpg" format from: 
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